Form-factors computation of Priedel oscillations in Luttinger liquids. 
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We show how to analytically determine for g < ^ the 
"Friedel oscillations" of charge density induced by a single 
impurity in a ID Luttinger liquid of spinless electrons. 

PACS numbers: 72.10.Fk, 73.40. Hm, 75.10.Fk. 



The general problem of a Luttinger liquid interacting 
with an impurity - that may have internal degrees of 
freedom - has attracted constant attention. A first rea- 
son is the wealth of physical applications: they include 
the anisotropic Kondo model, the double well problem 
and washboard potential problem of dissipative quan- 
tum mechanics, scattering through an impurity in quan- 
tum wires and tunneling through a point contact in 
the fractional quantum Hall effect Another reason 
is that this general problem is integrable, and therefore 
the possibility of obtaining exact solutions exists. Until 
recently however, such solutions had been restricted to 
thermodynamic properties. While of crucial experimen- 
tal interest, correlation functions and related transport 
and dynamical properties had remained inaccessible an- 
alytically. Numerical simulations were quite difficult, and 
not always conclusive. 

Recently, major progress has been made. Based on a 
new basis of massless quasiparticles suggested by inte- 
grability, together with a generalization of the Landauer 
Biittiker approach, DC properties have been exactly com- 
puted, in remarkable agreement with experimental re- 
sults 0]. Using the form-factors approach the 
dynamical properties of currents have also been obtained 
exactly 0|, or, more precisely, in closed forms that have 
an arbitrary accuracy all the way from the UV to the IR 
fixed points . 

The method used in 0| worked only for currents, ie 
for operators with no anomalous dimension. This was 
a major drawback, since many physical properties are 
described by more complicated operators. We show in 
this letter how operators with non trivial dimension can 
be also handled after some appropriate regularization. As 
an application, we determine the 2kp part of the charge 
density profile in a one dimensional Luttinger liquid away 
from an impurity, a problem which has attracted a lot of 
interest recently §, § 0. 

We start with the bosonised form of the model. The 
hamiltonian takes the form : 



H = 



dx [STT^tf -I- J_(a^^)2] + A COS 0(0), (1) 



where we have set vp = g. Then for the Friedel oscilla- 
tions, the charge density operator is just : 



p{x) = po + 2dx(j) H cos[lkpx + <t>{x)]. 



(2) 



with Po = ^ the background charge. We decompose 
this system into even and odd basis by decomposing 
4> = 4>L + 4>R and setting : 



1 



(fi" {x -\- t) = ^ [(f)Lix,t) + (l)B,{-X,t)] 
^°{X + t)^^ [(j)L{x, t) - M-X, t)] 



(3) 



Observe that these two field are left movers. We now fold 
the system by setting : 

= V2ip%x + t), X <0 (f>% = y/^cfy^i-x -^t),x <0 
(j)l^^/2ip''{x + t), X < (f>% = -V^c/y^i-x + t),x < (4) 

and introduce new fields 4)'^'° — 0^'° + (/'^°. The density 
oscillations now read : 

{p{x) - Po 



Po 



(t)°(x) (t>'^{x) 
cos(2fci?a; + ?7f)(cos — - — )(cos — - — ) 



(5) 



with rjF the additional phase shift coming from the uni- 
tary transformation to eliminate the forward scattering 
term. 0° is the odd field with Dirichlet boundary condi- 
tions, at the origin 4>°{0) — leading to : 



1 

(cos — ) cx 

2 \ a; 



ff/2 



(6) 



and the 4>'^ part is computed with the hamiltonian : 



= 1 f dx [87rgn^2 ^ ^(9^0^)2] ^ ^^^^ 

2 J-oo ^T^a 2 



On general grounds, we expect the scaling form : 
(cos^>« (-) F(Axi"^), 



(7) 



(8) 



where F is a scaling function to be determined. Note 
that even the small x behaviour of this function was no 
known in general. 
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Our approach is based on the fact that both systems 
are integrable. By considering the free boson as a Hmit of 
the sine-Gordon model we describe it using a basis 
of quasi particle states, the quasiparticles being massless 
solitons/antisolitons and breathers, with factorized scat- 
tering. The boundary interaction is then described by a 
scattering matrix, which is elastic , . 

To compute the correlation functions , it is convenient 
to represent the boundary interaction through a bound- 
ary state |B) : 



\B) = exp 



E 

.ei,e2 



/oo jn 
-oo ^''^ 



(9) 



where the e^'s denote the type of particles (solitons/anti- 
solitons or breathers) and the superscript denotes 
whether they are left or right movers since they are mass- 
less particles. Here Ob is a scale related to A encoding the 
boundary interaction, A — > corresponds to Ob — > — oo 
and A^cxDto^s-^-cxD. Kf^^^i is related to the reflection 
matrices jlj] . As usual we have used rapidity variables 
to encode energy and momentum. For solitons and and 
antisolitons for instance, e = ±p = /^e^, ^ an arbitrary 
energy scale. 

The one point function of interest reads then 
(O|cos^|i?). To use (^), we need first the matrix ele- 
ments of the operator cos ^ in the quasiparticle basis: 
these follow easily from the massive sine- Gordon form- 
factors 1^. Unfortunately, as discussed briefly in j^], the 
resulting integrals are all IR divergent! This was not the 
case for the current operator, whose form factor has the 
naive engineering dimension of an energy, leading to con- 
vergent integrals. Some sort of regularization is needed, 
and the correlations of cos ^ with a boundary had re- 
mained so far inaccessible. Our purpose is to show how 
to cure this problem. 

To explain our strategy, we consider first the case g — 
1/2. Here, the friedcl oscillations arc simply |jll[ related 
to the spin one point function in an Ising model with 
boundary magnetic field. By using the same approach 
as the one outlined before, one finds the following form- 
factors expansion : 



^ 1 /"OO 

))-E;^/ n l^tanh 

„_n J-OO 

n,;<j (tanh 



n=0 



(10) 



The integrals are all divergent at low energies, when Oi — s- 
— oo and the integrand goes to a constant. Let us then 
introduce an IR cut-off ( we chose > Omin and set 
A = e^'"'" ) and take the log of the previous expressions (a 
similar method has been used in [|j to study the UV limit 
of massive correlators. See also , |16 , 0). Ordering 



this log by increasing number of integrations, one can 
show that each term diverges as In A. Moreover, since 
the divergence occurs at very low energy, where the tanh 
goes to unity, the amplitudes of these In A do not depend 
on Ob (for Ob 7^ —00), ie on the boundary coupling. It is 
then easy to get rid of the cut-off: we simply substract 
the log of the IR spin function, ie we substract the same 
formal expression with 9b = 00. The first two terms of 
the resulting expression read : 



In- 
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(11) 



where we have set fj, = 1, m = e^S Tb = e^^ oc A^^^^ 




1 00 



FIG. 1. Accuracy of the finite Tb over the IR part of the 
envelope of p(x) for g — 1/2. 

Clearly, the integrals are now convergent at low ener- 
gies, and we can send A to zero. Since the IR value of the 
one point function is easily determined by other means, 
< a{x) >/fl'(x x~^^^ jl^, we can now obtain < a{x) >Tb 
from (|l^). Hence the procedure involves a double regu- 
larization. Of course, there remains an infinity of terms 
to sum over. However, as in the case of current oper- 
ators, the convergence of the form- factors expansion is 
very quick, and the first few terms are sufficient to get 
excellent accuracy all the way from UV to IR. To illus- 
trate this more precisely, we recall that for (7=1/2 ( pT| ) 
can be resummed in closed form, giving rise to : 

{<^{x))tb 



Re 



{a{x)) 



^ = ^^/2^e^^-KoixTB). (12) 

IR V"" 
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By reexponentianing the two first terms in (|fl|) , one gets 
a ratio differing from ( [l^ ) by at most 1/100 for xTb G 
[0,00) (see figure 1). 

By reexponentiating the first three terms, acuracy is 
improved to more than 1/1000. Clearly, the form-factors 
approach thus provides analytical expressions that can 
be considered as exact for most reasonable purposes. 

It is fair to mention however that, at any given order 
in (11), the exponent controlling the a; — *■ behaviour 
is not exactly reproduced, as could be seen on a log-log 
plot. For instance, the first term is immediately found to 
produce a behaviour R{x) cx x^/'^, to be compared with 
the result Rexact{x) oc x^^^lnx. The comparison of the 
exact result (|l^) and of (^l|) show that the form factors 
expression has, term by term, the correct asymptotic ex- 
pansion ie the IR expansion in powers of jy^- Adding 
terms with more form-factors simply gives a more ac- 
curate determination of the coefficients. This is to be 
compared with the results of [Q for eg the frequency de- 
pendent conductance, where the form factors expression 
had the correct functional dependence both in the UV 
and in the IR. This is not to say that our method is in- 
efficient in the UV, because we know, at least formally, 
all the terms. In fact, we will show in what follows how 
( pT| ) can always be resummed in the UV, and that the 
exponent can be exactly obtained from our approach too. 

The regularization is the same for other values of g. 
Here, we discuss g — 1/t with t integer. For these values, 
the scattering is diagonal and the form factors are rather 
simple. To obtain them, we take the massless limit of the 
results in j|] and impose that half of the quasiparticles 
become right movers and half become left movers, since 
the boundary state always involve pairs of right and left 
moving particles. It is in fact easier to take that limit 
if we change basis from the solitons and anti-solitons to 
> ±1^ >). In that case, the boundary scattering 
matrix become diagonal and the isotopic indices always 
come in pairs. The reflection matrices in this new basis 
are given by : 

K^{e) = -e^f (2-*) tanh(^^ + i^^^)R{i^ - 9) 

K+ie) =e'^^^~'^Rii^-e) (13) 



with 



R{9) = exp L 



2y sinh(2y) cosh((t — 1)?;) 



(14) 



The breathers reflection matrices are given in |19|. 

The case g = 1/2 has already been worked out, so 
let us concentrate on g — 1/3 as an example. Then, in 
addition to the soliton and anti-soliton, there is also one 
breather. The first contribution to the one point function 
comes from the two breathers form- factor, with one right 



moving and one left moving breather. It is given by a 
constant : 



LR 



Cl, 



(15) 



and this obviously leads to IR divergences. Other con- 
tributions come from 2n breathers form-factors, and 4n 
solitons form-factors The whole expression can be con- 
trolled as for g — 1/2, by taking the log, and factoring 
out the IR part. Setting c{x) = cos we organize the 
sum as follows : 



In 



i44r^=lni?(^)+lni?(^) + ... 



— -iii.f - -riii.f - -r--- (16) 

with the subscript denoting the number of intermediate 
excitations. 

Then, using the explicit expressions for g = 1/3 we 
find : 



lni?(2) ^ 2cie2^^«^£;i(-2\/2Tsa;), 



(17) 



where Ei is the standard exponential integral. The next 
term Ini?*^^^ is a bit bulky to be written here, but it 
is very easy to obtain - similar expressions have been 
explicitely given in |Q.This is all what is needed for an 
accuracy better than 1 percent. In figure 1 we present 
the results of the ratio at g = 1/2, 1/3, 1/4 for the Friedel 
oscillations. It should be noted that this ratio is just the 
pinning function of reference Q and our results agree well 
qualitatively with the results found there. 

As mentioned before, the deep UV behaviour is a little 
more difficult to obtain: the accuracy is good because the 
ratio goes to zero anyway, but the numerical evaluation of 
the power law would not be too accurate with the number 
of terms we consider. Fortunately, the full form-factors 
expansion allows the analytic determination of this expo- 
nent. First, observe for instance that in (|l^) the integrals 
converge for all Tb ^ 0, but strictly at Tb = 0, they do 
not. To get the dependence of (c(x)) as Tb — > 0, we 

Ic(x))t 

will consider, the logarithm of another ratio. In , ) ,\\ ^ , 

where x and x' are two arbitrary coordinates. For this 
ratio, even at Tb = 0, the integrals are convergent. But 
Tb = is the UV fixed point, with Neumann boundary 
conditions. While the one point function {c{x))uv van- 
ishes, the ratio of two such one point functions is well 
defined, and can be computed by putting an IR cut-off 
(a finite system). One finds that it goes as [x/x'y/'^. By 
regularity as Tb ^0, the same is true for the ratio close 
to Tb = 0, and thus one has 



{c{x)) o^{xTBy/^x{TB)-^0 



(18) 



This shows that the universal scaling function in (jg) be- 
haves as F{y) cx for g < ^. This exponent can 
actually be obtained by perturbation theory. Indeed, the 
first term in the perturbative expansion of (c(x)) is 



3 



dy 



(x2 + y2)3 



(19) 



For g < 1/2, this integral diverges in the IR. To regu- 
late it, we need to put a new cut-ofF: since there is no 
other length scale in the problem, this can be nothing 
but I/Tb- Changing variables, the leading behaviour is 
x^^^Tg (X a;^/^A^/^~^, in agreement with the previous 
discussion. The exponent coincides with the result of the 
self-consistent harmonic approximation |8|; but it is im- 
portant to stress that the latter is vahd only foY g « 1. 
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FIG. 2. Ratio of the finite Tb over the IR part of the 
envelope of p{x). 

The function F{y) behaves aylny for 5 = 1/2. For g > 
1/2, its behaviour is simply F(y) cx y, as can be easily 
shown since the perturbative approach is now convergent. 
As we approach g = 1/2 this exponent seems to become 
asymptotic and is more difficult to get numerically ||]. 

The method presented here is very successful to get 
analytical results for g < 1/2 - although we limited our- 
selves to g — 1/t with t integer, all values of 17 < 1/2 are 
accessible, but computations are more complicated since 
the bulk scattering is non diagonal. The method should 
be generalizable to other problems, in particular the de- 
termination of the screening cloud in the anisotropic 
Kondo model , as will be reported elsewhere. The re- 
gion f; > 1/2 presents additional difficulties, unresolved 
for the moment - in particular, the massless limit of the 
form factors does not seem to be meaningful. Of course 
the case g — 1 can be solved by fermionization |^ . In our 
approach, this point is non trivial because of the folding. 
This folding however is necessary for any value g I: 
except at (7 = 1, the problem on the whole line would not 
be integrable. 
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